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ABSTRACT 


Shock tube experiments have been performed to 
de t e rm i ne the response of a hot-film sensor, mounted f 1 u s h 
on the side-wall of a shock tube, to unsteady flow behind a 
normal shock wave. The present experiments attempt to 
isolate the response of the anemometer due only to the 
change in convective heat transfer at the hot-film surface. 
The experiments, performed at low supersonic shock speeds 
in air, are described along with the data acquisition 
procedure. The change in convective heat transfer is 
deduced from the data and the results are compared with 
those from transient boundary- layer theory and another set 
of experimental results. Finally, a transient local heat 
transfer coefficient is formulated for use as the forcing 
function in a hot-film sensor instrument model simulation. 
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CHAPTER 1 
INTRODUCTION 

The transient response of a constant temperature hot- 
film sensor, mounted on the side wall of a shock tube, due 
to the passage of a moving normal shock wave and the 
ensuing convective heat transfer is analyzed. Also 

included is an interpretation of transient, compressible 
boundary- layer theory for application to shock tube 
experiments with heat convected from a side wall mounted 
hot-film sensor. This type of sensor is currently in use 
by many researchers doing experimental wind tunnel 
research. The effective use of this sensor requires a 
knowledge of its steady-state, as well as its transient 
response characteristics due to changes in flowfield 
conditions. The experimental steady-state response 

characteristics have been documented by Wusk , et aj. [1] ■ 

To obtain the transient response characteristics, a shock 
tube is used to induce an unsteady, compressible boundary- 
layer "flow across “the hot-'film sensor. 


'Numbers in brackets refer to the list of references. 
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1 . 1 Mot i vat i on 

Hot-film sensors are currently be i ng used 
experimentally for many purposes, including the detection 
of transition from laminar to turbulent flow in a flowfield 
[2]. The research presented in this thesis was originally 
motivated by a desire to use a hot-film anemometer in the 
detection and measurement of cross-flow vortices [l]. 1 11 

order to solve this problem, an understanding of' the 
transient response characteristics of the hot-film sensor- 
must be known. Consequently, this work is limited to 
analyzing the transient response of the hot-film sensor due 
to an almost instantaneous step-change in convective flow 
conditions. Furthermore, the research is directed by the 
desire to correlate the experimental results with existing 
transient boundary- 1 aye r theory [3] . This work will be oi 
interest to those using hot-film sensors for various 
applications, because a knowledge of the sensor’s steady 
and transient response assures a more accurate means of 
characterizing flow over a surface. 

1.2 General Problem Description 
The general problem treated in this thesis is to 
analyze the transient response of a constant-temperature 
hot-film sensor from shock tube experiments in order to 
determine how the convective heat transfer varies with 
time, and correlate these results with existing transient 
boundary- 1 aye r theory. Consequently, the hot-film sensor 
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response can then be predicted under variable flow 

cond i t i one . 

Normal shock theory indicates that a step change in 
pressure, temperature, and velocity occurs across the wave. 

A sudden change in flow properties causes a transient 
response in the hot-film sensor. Because the normal shock 
wave is moving in the laboratory frame of reference, a 
transient boundary layer develops, and the hot-film sensor 
responds to the change in convective heat transfer 

occurring at the sensor with a change in the anemometer 
output voltage. The subsequent boundary- 1 aye r build-up 

after the passage of the initial normal shock wave 

indicates that the boundary- 1 aye r thickness increases, 
which in turn decreases the convective transport of energy 

from the sensor as time increases. Theoretical work in the 

determination of the compressible boundary layer behind a 
moving normal shock has been done by Mirels [3,4] . 

1.3 Relevant Literature 

The development of the theoretical prediction of the 
transient convective heat transfer behind a moving normal 
shock was accomplished in two papers published by Mirels 

[3,4] . The first paper predicts the laminar velocity 
boundary- 1 aye r characteristics as well as the thermal 
boundary- layer characteristics behind a moving normal shock 
advancing into a stationary fluid [4] . The second of thes< 
two papers predicts the thermal and velocity boundary layer 
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profiles for both laminar and turbulent, flow behind a shock 
or thin expansion wave advancing into a stationary f 1 u i r \ 
[3]. These works are used extensively in this thesis as 
the basis for predicting the theoretical velocity and 
thermal boundary layer profiles across the hot -film 

sensor. DISA, a Denmark manufacturer of thermal anemomet ry 
instrumentation, produced results for hot-wire anemometers 
which qualitatively corroborate with the anemometer 
response for the shock tube tests done in this thesis 
research [[5] . Wusk , Carraway , and Holmes use an arrayed 
hot-film sensor for laminar bou ndary- 1 aye r studies aimed 
toward the detection of cross-flow vortices [1 ] . 
Furthermore, the same type of constant temperature hot-film 
sensor, which is employed in the present research, was 
tested by Wusk, et a ! . The arrayed hot - f i lm sensor' was 
mounted on a NASA NLF(1)-0414 nat u ral 1 am inar f 1 ow a i rfoi 1 
and tested in the NASA Langley Instrument Research Division 
(IRD) small calibration facility. The goal of that 

research was to take the initial steps in the steady-state? 
calibration of the hot-film sensor under constant flow 
conditions for the detection of spanwisc variations of heat 
transfer. In the DISA work, a hot-wire anemometer was 
moved radially from the shock tube axis to within 0.05 mm 
°f the wal 1 of the shock tube to measure the shock f rout, 
curvature. Also, a comparison of the experimental results 
from this thesis is made with the work by Davies and 
Bernstein [6] . Tn the experimental investigation by 
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Davies and Bernstein, a semi-infinite flat plate was 
mounted in the low-pressure chamber of a shock tube in 
order to determine the convective heat transfer rate due to 
a shock- induced boundary layer. Furthermore, Davies and 
Bernstein also use M i re 1 s ’ works [3,4] for correlation with 
their experimental measurements. 

A paper published by Roberts, et aj_ [13] used high 

frequency response (typically of order 10 6 Hz) miniature 
surface thermocouples mounted on the surface of the 
Austrailian National University (ANU) T3 shock tube to 
detect convective heat transfer and compare the heat 
transfer on ’’clean” and ’’dusty” surfaces. The research had 
test times on the order of 200 /is due to the very high 

pressures at which the tests were conducted. Roberts, et 
aj_ illustrates a comparison of his results with that of 
Mirels’ turbulent boundary- 1 aye r theory [3]. 

1.4 Problem Statement 

The problem to be analyzed in this thesis is the 

variation of convective heat transfer with time due to 
thermal and velocity boundary- 1 aye r growth behind a moving 
normal shock wave across a constant-temperature hot-film 
sensor. A shock tube is used as the experimental means of 
analyzing this problem. As the shock passes over the 
sensor, a flow is induced behind the normal shock wave 
which causes the convective heat transfer above the hot- 
film sensor to change, and thus produces a change in the 



anemometer output voltage. The purpose of this thesis is 
to measure the transient convective heat transfer using the 
surface temperature and heat transfer rate provided by the 
hot-film sensor. Results are compared with experimental 
results of other researchers, and correlated with 
transient, compressible bou ndary- 1 aye r theory. 

In the following chapter, the transient, compressible 
boundary- 1 aye r theory is presented and briefly discussed. 
The theoretical formulations applied in this thesis are 
interpretations taken from M i re 1 s * works [3,4j regarding 
the prediction of the boundary- 1 aye r growth behind a moving 
normal shock wave. The following chapter contains the 
reference frame transformation used to modify the governing 
equations into a more useful form. Also presented are the 
laminar and turbulent correlations for the* convective heat 


transfer . 



CHAPTER 2 

THEORETICAL FORMULATIONS 


2.1 Reference Frame Transformations 

A boundary layer is established along the shock tube 
wall and across the hot-film sensor as the normal shock 
wave passes. The boundary- 1 aye r growth affects the 

convective heat transfer rate from the sensor. The normal 
shock wave is assumed to travel at a. constant velocity, 
parallel to the wall, into a stationary fluid. The choice 
of coordinate systems for the theoretical analysis is 

considered here in some detail . 

A coordinate system, x and y, is fixed with respect 
to the shock tube wall. The corresponding velocities are 
u, v. In this coordinate system, the wave is moving at a 
constant velocity and is considered to move in a time 
dependent (i.e. unsteady) reference frame. In order to 

make simplifications in the governing equations, the 

application of a steady reference frame is required. A 
coordinate system is employed in which the wave is 

considered stationary and the wall is moving at a constant 
velocity equal in magnitude to the shock wave velocity in 
the unsteady reference frame. In this reference frame, x 
and y represent the coordinates fixed to and moving with 


7 



8 


the wave. The corresponding velocities in t li i s reference 
frame are u and v. By selecting this reference frame , the 
coordinate system is moving with the constant wave velocity 
and the flow is considered steady. Figure 2.1 illustrates 
the coordinate systems used in the theoretical analysis. 

The following assumptions are employed in order to 
develop the governing equations in a more useful form. 
Because the coordinate system is moving with the wave , the 
time derivative terms in the governing equations are zero. 
With laminar bou ndary- 1 aye r flow, retaining compressibility 
effects and assuming dP/dx = 0, the Prandtl boundary- 1 aye r 
equations apply for flow near the wall [7]. By assuming 
the wave to travel at a constant velocity, the velocity of 
the flow behind the wave is also traveling at a constant 
velocity. Because of this assumption and the fact that the 
sensor is treated as a flat plate, Bernoulli’s equation is 
employed to validate the assumption dP/dx = 0. 
Consequently, the governing equations are: 

Cont i nu i ty , 


d(pu) 3(pv) 

dx dy 


( 2 . 1 ) 


Momentum , 


u |y + 

ax ay 


1 d_ 

P dy 


0u1 


P 3* 


(2.2) 


Energy , 


pc p 


r„OT , w dTi 

_ d 

f, dT^ 



dy 

L k spj 

+ "I^J 


(2.3) 
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- ig-re 2 i. Coora -ate Systems Used in 
Boundary Layer Anaiys.s 



Idea.! Gas Law , 


1 () 


p — pRT - const. (‘2.4 ) 

subject to: (x > 0) 

u(x , 0) — U W u(x , oc) = u e 

v(x,0) = 0 (2.5) 

T(x , 0) = T w T(x , oo) = T e 

2.2 Similarity Transformation of Governing Equations 

The governing equations can be transformed to a form 
that is easily integrated numerically. From Mi rel s [3] , 
variable transport and thermal properties are applied in 
which /i and k are assumed proportional to T. Also, Cp and 
<7 are assumed to be independent of T, but are evaluated at 
the wall temperature, T w . The thermal and transport 

properties are ar bitrari ly referenced to the wa) 1 . 
Defining the similarity parameter as: 


V 


N 


u e 

2xi^ w 


J 

I 

0 


T(x,y) 


dy 


( 2 . 0 ) 


the governing equations can be transformed by the usual 
boundary layer similarity transformations. Thus, the non- 
linear momentum equation becomes: 

f + ff" =0 (2.7) 


subject to: 


f(0) = 0 
f'(0) = ^ 
f'(oo) = 


( 2 . 8 ) 


1 



linear and can be 


The transformed energy equation is 
expressed as a linear superposition of the solutions for 
the case with zero heat transfer plus the case with heat 
transfer [3]. Consequently, the energy equation becomes: 


T t T u w if u e r (*») r 

TI = 1 + N" 1 J2T7c p , w f 


T w _ It 

TV T e 


s(»?) 


(2.9) 


where the functions r(,) and s(r,) are determined as follows: 
Zero Heat Transfer Case, 


r" + Mr' = ( f 'V 


Ps- 1 ] 


( 2 . 10 ) 


subject, to : 


'(oo) = r'(0) = o 


and Heat Transfer Case 


i n + u w f = 0 


( 2 . 11 ) 

( 2 . 12 ) 


subject to : 

S(0) = 1 (2.13) 

e(oo) = 0 

High speed convective heat transfer is accomplished 
respect to a recovery temperature, T r , which takes the 

f o rm : 

i . r u w ,f u j r (°) 

1 + 2T e c p ,w 


T r _ 


(2.14) 


With the energy equation defined by equation (2.9), 
further reduction may be expressed in terms of the 
convective heat flux at the shock tube wall: 



n 

Qw 


(Tw-T r ) 


1 


k dj 

dy 


kwB,(0 \f5§ 


(2 . 15,i 


This equation becomes the most useful form in comparing the 
experimental data to the theory for laminar f 1 ow 
cond i t i ons . 


2.3 Laminar Correlation 

The laminar correlation for the convective heat 
transfer occurring above the hot-film sensor is deve 1 oped 
from the governing equations using the similarity and 
reference frame transformations. Because radiative and 

conductive heat transfer effects are neglected, the 
resulting equation for the heat transfer from the hot-film 
sensor is that of convective heat transfer. By definition, 
the convective heat flux in high velocity boundary- layer 
flow is defined: 

q" = h(T w -T r ) (2.16) 


A direct substitution with Eq . (2.15) results in an 

expression for the convective heat transfer coefficient, h. 
Due to the reference frame transformation, Eq . (2.15) is 

val id for a coordinate system moving with the shock wave 
velocity. Since the experimental results are gathered in 
the laboratory (i.e. unsteady) frame of reference, the 
theoretical equations are transformed to the lab reference 
frame to maintain consistency between theoretical and 
experimental resu 1 *^. Consequently, the equations are 



13 


initially derived in the wave reference frame and 
transformed to the lab reference frame by initiating the 
following substitutions for position and velocity, 
respect i vely : 


x = u w t — x 
u e — U W 

By utilizing these relations, the convective 
coeff ic ient becomes : 


( 2 . 17 ) 

( 2 . 18 ) 
heat transfer 


h 



(2 . 19) 


The convective heat transfer coefficient becomes a function 
only of time, t. 

As stated in Section 1.4, an objective of this thesis 
is to compare the results of this thesis to the work of 
Davies and Bernstein [6]. In order to make this 
comparison, the Reynolds, Nusselt, and Stanton numbers are 
defined. M i re 1 s [3] defines a Reynolds number with respect 
to a fluid particle traveling behind an expansion wave in a 
stationary (lab) reference fraine. Applying Mirels’ 
definition of Reynolds number to a shock wave and employing 
Eqs. (2.17) and (2.18), a local Reynolds number for the 
fluid particle behind the wave in the lab frame of 
reference is defined: 


R- e xp 


( u w 


-u e ) 2 t 

‘'w 


( 2 . 20 ) 


Davies and Bernstein [6] illustrate their results using a 


Reynolds number defined wit,h respect to a fluid particle 
traveling with the shock velocity. The fol lowing 
substitution transforms Eq . (2.20) into a Reynolds number 
defined with respect to a fluid particle traveling with the 
shock velocity: 


R e xs — ^ e xp^tT^] 


( 2 . 21 ) 


The Nusselt number is a non-dimensional temperature 
gradient, when surface and free stream temperatures are 
fixed, and provides a measure of the convective heat 
transfer occurring at the surface. Since a change in 
reference frame is used to transform the governing 
equations into steady boundary- 1 aye r equations, the Nusselt 
and Stanton numbers are therefore defined in this steady 
reference frame and then transformed into the laboratory 
(unsteady) reference frame. By definition, the Nusselt. 
number, with respect to the shock wave reference frame, is 
defined as: 

Nu = ^ ( 2 . 22) 

K w 


After substituting Eqs . (2.17) and (2.19) into Eq . (2.22), 

the Nusselt number takes the form: 


As an alternative 
dimensionless conv 



to using the Nusselt 
ective heat transfer 


(2.23) 

number for the 
parameter, the 


i s 


Stanton number (modified Nusselt number) 


used . 


The 
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derived Stanton number takes the form: 


St = 


-s'(0) 


/>|2 1 jRe xp 


(2.24) 


This equation is then applied to compare this work with the 
work of Davies and Bernstein [6]. These derivations are 
necessary for laminar flow conditions; however, M i re 1 s [3] 
develops a turbulent correlation that is presented in the 
fol lowi ng sect, ion . 

2.4 Turbulent Correlation 

The turbulent boundary- 1 aye r solution for convective 
heat transfer behind a moving expansion wave has been 
developed in detail by Mirels [3]. The solution is 

obtained for a moving shock wave by extending empirical, 
semi - inf inite flat plate boundary-layer theory to the case 
of a moving wall. Mirels [3] assumes a (l/7) th power 
velocity profile and extends this profile to a moving wall. 
Also, compressible turbulent flow over a sem i - i nf i n i te flat 
plate is approximated by evaluating the fluid properties at 
a mean static temperature, T m , for the Blasius relation of 
incompressible turbulent flow past a sem i - i nf i n i te flat 
plate [3]. A reasonable estimate for the mean temperature, 
T m is given by Eckert [11] : 


T m = 0 . 5 (T w +T e ) + 0.22(T r — T c ) 


(2.25) 


The final expression of the integral form of the momentum 
equation, after appling the Blasius relation, takes the 



■form : 


1 (i 


r w 

Pe u e 


0 . O46O-(0/6 ) 


4V 


« j 


( 6 / 6 ) 


Uw 

Ue' 


3 

.isr 




( 2 . 20 ; 


whe re , 



(2 .27 j 


and (0/£) is the ratio of the boundary-layer moment urn 

thickness to the fluid velocity boundary- 1 aye r thickness 

[3, defined by Eq . (46a)] . 

The Reynolds-Col burn analogy [12]: 

St<r3 - ^ (2.28) 


is applied with the skin friction coefficient defined with 
respect to the lab reference frame as: 

c.c = t — (2.29) 


f |pw[u w -u e ] 5 


and substituting Eqs. (2.26) , (2.28) and (2.29) , the final 

equation for the convective heat transfer behind a moving 

normal shock in a compressible , turbulent, boundary- 1 aye r 
flow takes the form: 


St Re 


i 

xps 


= 0 . 0460* 



^m 


2 

3 


(2.30) 


Thus, a comparison of the theory for a compressible, 
unsteady, turbulent-boundary layer flow can be made with 
the work by Davies and Bernstein [6] as well as Roberts, et 
al [13] , and with experimental results developed here. 




and applying Eqs . (2.20) , (2.22) , and (2.30) , a turbulent 
heat transfer coefficient is derived. Furthermore, M i re 1 s 
[3] neglects the variation of fluid state properties which 
results in a simplification of Eq . (2.30). The turbulent 
convective heat transfer coefficient takes the form. 


h 


i 3 

0 . 0460* k m <7ni ( u w “ u e) 
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* 

6 



(2.32) 


The experimental procedures and set-up for the shock 
tube tests conducted using the NASA Langley Shock Tube are 
discussed in the following chapter. Also presented are the 
equations used to reduce the experimental data into a form 
for comparison with the theoretical results and with 
experimental work done by Davies and Bernstein [6] and 
Roberts , et aJL [13] . 



CHAPTER 3 

SHOCK TUBE EXPERIMENT 


3 . 1 Pu rpose 

The shock tube is utilized as the test apparatus for 
the analysis of the transient response of a hot-film sensor 
on a foam substrate mounted flush in the shock tube wal 1 . 
The shock tube is utilized as the test apparatus because an 
almost instantaneous step-change in flow properties is 
obtained across the moving normal shock as the shock passes 
over the sensor, and thereby causes a transient boundary- 
layer build-up. Also, the unconstrained access and 
simplicity in operating the shock tube cannot be overlooked 
as a driving force in the decision to utilize the shock 
tube as the test apparatus for transient response 
experiments. The ultimate goal behind the experiment is to 
determine the hot-film anemometer response due to 
variations in convective flow conditions and to ascertain 
the repeatability of the experimental results. 

3.2 Experimental Set-up 

The equipment utilized in the test is listed in 
Appendix A. The actual shock tube assembly used for the 
tests is seen in Fig. 3.1. The experimental equipment 
configuration outlining the flow of experimental 
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I’i^ure 3.1 . Shock Tube Assembly 
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information is seen schematically in Fig. 3.2. In 
referring to Fig. 3.2, note that the shock tube contains 
three sets of downstream ports on the tube wall . Two PCB 
112A21 Piezotronic H i gh - Reso 1 u t i on Pressure Transducers are 
placed in the top ports at positions 1 and 2 of the shock 
tube. A pressure transducer is mounted in the top port at. 
position 2 to give an accurate measure of the speed of the 
shock wave as the wave travels down the shock tube. By 
noting the time-synched transducer responses displayed by 
the Gould Digital Oscilloscope, measurements are obtained 
for the time interval during which the shock wave travels 
from position 2 to position 1 of the shock tube. By 
knowing the distance between positions 1 and 2 (24.0 in or 
61.0 cm), the actual speed of the shock wave is calculated. 
Consequently, this calculation gives a means to verify the 
predicted theoretical wave speed determined from normal 
shock equations [8] . 

The hot-film sensor is placed on an insulating foam 
substrate and mounted on a plug placed in the position 1 
side wall port of the shock tube. The physical dimensions 
of the sensor are illustrated in Fig. 3.3. Judge [9] gives 
a listing of the material properties of the hot-film sensor* 
and substrate. The anemometer and pressure transducer 
signals at position 1 of the shock tube are synchronized in 
time in order to show the simultaneous variation in the 
anemometer response and pressure transducer response. A 
Type T thermocouple is inserted into the top port at 
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position 0 of the shock tube. The thermocouple oxtonds 0.5 
in (1.3 cm) from the wall into the flow and measures the 
ambient air temperature in the shock tube downstream of the 
anemometer near the endplate (Fig. 3.2). Also, a Type T 
thermocouple is mounted on the foam substrate beside the 
hot-film sensor. Figure 3.4 shows both the hot-film sensor 
and the thermocouple mounted on the foam substrate. The 
thermocouple measures the variation in the substrate 
temperature due to the constant elevated operating 

temperature at which the hot-film sensor is maintained. 
Initially, the hot-film sensor is maintained at an elevated 
temperature of 380 K which is 85 K above ambient. This 
corresponds to an overheat ratio of 1.3. The response of 

the anemometer is also tested for overheat ratios of 1.4 
and 1.5 (see Section 3.3 for anemometer temperature 
calibrations). The overheat ratio is defined as the ratio 
of the heated sensor resistance to the cold sensor 

resistance. The effect of increasing the overheat ratio is 
to increase the anemometer sensitivity to mass flow 
fluctuations and increase the frequency response of the 
anemometer. The substrate temperature, measured at the 

thermocouple location, maintains approximately the same 
average temperature (22.8 ± 0.6 ‘C) throughout the series 

of tests for all three overheat ratios. 

The pressure transducers are powered by two PCB 484B 
Line Power Units and the responses are displayed on a Gould 
Digital Storage Oscilloscope, Type 4035. The pressure 
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t, ransd u ce r at, position 1 and the anonuiinclcr are conru ct ( d 
to channels 1 and 2, respectively of a Tektronix Digital 
Oscilloscope, Model 2430A . A personal computer is used as 
a controller for this experiment in order to initialize the 
Tektronix Oscilloscope parameters necessary in recording 
the anemometer and channel 1 pressure transducer data. A 
program written in GW-BASIC is used to communicate between 
the computer and oscilloscope via an IEEE-488 Interface. 
Appendix B contains a listing of the data acquisition 
program used in this experiment. For each test, the 

oscilloscope stores the voltage values of the anemometer 
response, which correspond to changes in the convective 
heat transfer across the sensor, and transfers this 
information onto a floppy disk for permanent record and 


future data analysis. The pressure transducer response at 
position 1 is also recorded onto a floppy disk. The 

pressure transducer and hot-film sensor responses recorded 
and displayed by each of the digital oscilloscopes are 
plotted for a hard copy of the output. The anemometer and 
channel 1 pressure transducer responses are plotted on a 
Hewlett-Packard 7470A Digital X-Y Plotter and both pressure 
transducer responses are plotted on a Hewlett-Packard 7046A 
X-Y Recorder. Also, a barometer is placed in the 

laboratory to obtain a more accurate measure of the ambient 


pressure 


at- the time o*f testing. 
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3 . 3 Kxpe r i men ta L Procod u re 

The shock tube tests are completed in a series of 
steps for a range of compression chamber (driver) 
pressures. The overheat ratio of the sensor is varied from 
1-3 to 1.5 for each range of driver pressures. Also tested 
is the effect of the hot- film sensor response to changes in 
sensor orientation at the same driver pressure. 

The initial set-up parameters for each test are 
recorded prior to the breakage of the Mylar diaphragm. 
These parameters are: (1) ambient temperature in the shock 
tube, (2) atmospheric pressure, and (3) substrate 
temperature. The computer adjusts the scope settings 
before each test and is used to record the fluctuating 
output voltage of the hot-film sensor due to the changes in 
convective heat transfer occurring across the sensor. 

Before the initial testing, the anemometer’s frequency 
response is tested using a 30 kHz square wave signal . A 
signal response of approximately 93 kHz is obtained with a 
sensor overheat ratio of 1 .3. The overheat, ratio of 1 .3 
corresponds to a total heated sensor resistance of 14.62 ft 
computed by taking the sum of the shorted sensor resistance 
at ambient temperature multiplied by the overheat ratio 
plus the cable resistance. The cable resistance is 
frequently calibrated during the experiment and is found to 
be very stable. Also, the sensor - reliably maintains its 
temperatu re - res i stance calibration, and consequently its 

i * -quency response for the tests at each 


sensitivity and 
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particular overheat ratio. At an overheat ratio of 1.4, 
the total heated sensor resistance is 15.56 U with a 
frequency response of approximately 104 kHz. Also, the 

frequency response is approximately 105 kHz with a total 
heated sensor resistance of 16.56 0 for an overheat ratio 

of 1.5. 

The hot-film sensor is postcal i brated by WYLE 
Laboratories in conjunction with NASA Langley IRD. The 

calibration produces an indication of the temperature 
the hot-film sensor at each overheat ratio. The 

calibration is performed by placing the sensor, mounted on 
the insulating foam substrate plug, into a temperature 
controlled air oven. The oven is elevated through a range 
of temperatures from approximately 15 *C to 150 C and 
subsequently decreased in temperature through the same 
range. Using a 4-wire resistance measurement, the 

resistance of each sensor is recorded corresponding to the 
steady state temperature of the oven for each data point. 
The plug consists of four sensors mounted on the foam 
substrate (see Fig. 3.4). Sensor #4 is shorted and is used 
to give an indication of the lead resistance for the 
sensors. Only sensor #1 is used throughout the series ol 
shock tube tests. Once the calibration is completed, a 
second order polynomial fit is placed through the data. 
Thus, an equation for the hot-film sensor calibration is 
obtained in which the hot-film sensor operating temperature 
is a function of the sensor resistance. Figure 3.5 



illustrat.es the obtained calibration curve; and, for ri 
specif ied sensor resistance at the time of testing for the 
corresponding overheat ratio, the operating temperature o f 
the hot-film sensor can be obtained. The non-linear curves 
t h rou gh the data are least squares best fits. 

Before each test, the end plate of the shock tube is 
removed and the debris remaining inside the tube due to t h< 
diaphragm rupture is blown out. The shock tube is cleaned 
after each test in order to minimize the amount of flow 
disturbance present in the shock tube. A Mylar diaphragm 
is inserted into its holder and placed in position in the 
shock tube (see Fig. 3.2). The Mylar diaphragm is cut and 
shaped from a sheet of Mylar with thickness corresponding 
to an estimated driver rupture pressure. As the driver 
rupture pressure is increased, the thickness of the Mylar 
diaphragm is also increased. The computer program is 
implemented to initialize the Tektronix Oscilloscope 
parameters and record the anemometer and position 1 
pressure transducer response data. After these procedures 
are completed, the driver section is pressurized until the 
diaphragm ruptures. At the time of the diaphragm rupture, 
the maximum driver pressure attained from a Wallace and 
Tierman pressure gauge is recorded, and the oscilloscope 
captures the transient anemometer and position 1 pressure 
transducer response data as the shock passes the hot-f i lm 
sensor . 

In order tc capture the pressure transducer response 
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Figure 3.5. Hot-Film Sensor Calibration 
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and hot-film sensor response, the digital oscilloscopes are 
set in a pretrigger mode. The pretrigger is used to 
capture the initial signal of the flow disturbance prior to 
the passage of the shock wave. Consequently, the entire 
response of the anemometer and pressure transducer can be 
displayed and recorded by the oscilloscope. The Tektronix 
Oscilloscope, which captures the anemometer and position 1 
pressure transducer signal, is initialized in the single 
sequence mode with a trigger level on channel 1 of 100 mV. 
As the shock passes position 1, the pressure increase due 
to the physical discontinuity of flow properties across the 
wave causes the scope to trigger. Once the 100 mV increase 
in voltage is obtained, which is almost instantaneously, 
the oscilloscope records the transient response of the 
pressure transducer as well as the anemometer. The Gould 
Oscilloscope, which records both pressure transducer 
responses, accomplishes the same task but the trigger level 
is adjusted manually to trigger off the positive rising 
edge of the pressure pulse. 

The tests are conducted at 4 driver gauge pressures 
of approximately: (1) 26 psig (179 kPa) , (2) 49 psig (338 
kPa) , (3) 56 psig (386 kPa) , to (4) 79 psig (545 kPa) . At 
each of these driver pressures, the test is repeated for 
which the time base of the oscilloscope is changed from 500 
/is/div to 20 ps/div. The sequence of tests is performed at 
an overheat ratio of 1.3 and repeated for overheat ratios 
of 1.4 and 1.5. Once these tests are completed, the 
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effect, of the orientation of the sensor is tested. The 
sensor is rotated at various angles relative to the 
standard operating position at a driver pressure of 
approximately 26 ps i g (179 kPa) and at an overheat ratio of 
1 . 3 (see Fig. 3.-4 for standard ope rat i ng pos i t ion) . 

3.4 Data Reduction 

In order to make a comparison of the experimental 
results with theoretical results, the experimental data is 
converted into a more useful form. First note that the 
effects of heat conduction into the substrate are neglected 
since the frequency response of the substrate will not 
allow the substrate to respond during the approximate 3.0 
ms test time. It is assumed that the hot-film sensor has 
been operating long enough to reach a steady-state 
condition. For the tests, the anemometer is switched 
operational approximately 15-30 minutes before the first 
test; long enough for the substrate heating to stabilize. 
The stabilization of the substrate heating is observed by 
monitoring the substrate thermocouple. It is also assumed 
that radiative heat transfer effects from the sensor are 
steady during the shock passage. Consequently, only the 
change in heat transfer from the hot-film sensor due to 
convection is observed during experiment times which are 
less than 3.0 ms. 

To fully understand the development in transforming 
the experimental data into a more useful form, it is 



necessary to understand the concept for the operation of 
hot-film sensors. The principle behind constant 

temperature anemometers is relatively simple. The 

electronic feedback system of the anemometer maintains the 
sensor temperature at a constant value . As a moving fluid 
convectively cools the hot-film sensor , the electronic 
feedback system of the anemometer increases the current 
supplied to the sensor in order to maintain a constant 
temperature. The sensor is actually one leg of a 

Wheatstone bridge. The sensor’s electrical resistance is 
proportional to temperature, and the resistance becomes a 
measure of sensor temperature [10] . A resistance change of 
the sensor is compensated by the control circuit sending 
enough current (Is) through the hot-film sensor to restore 
the bridge to its original set value. Consequently, the 
change in current to the sensor (Is) determines the 


electrical 

powe r 

d i ss i pated 

by 

the 

senso r , 

and 

the power 

increase be 

comes 

a measure 

of 

the 

rate of 

heat transfer 

from the 

sensor 

Due to 

the 

assumpt ions 

of 

negl igi bl e 


changes in conduction and radiation during short experiment 
times, the power loss is a measure of the convective heat 
transfer from the hot-film sensor. 

The experimental data is transformed into a more 
useful form by relating the output voltage of the hot-film 
sensor to the convective heat transfer. A series of steps 
is followed in order to transform the data. In order to 


determine the sensv r current, both the mean voltage output. 
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measured at each overheat, ratio, and the fluctuating output 
voltage due to the hot-film sensor reponso are used. At an 
overheat of 1.3, the output voltage level is measured to be 
2.85 V. Also, at overheat ratios of 1.4 and 1.5, the 
output voltage values are measured as 3.17 V and 3.43 V, 
respectively. The following equation is used to determine 


the sensor current,: 



(3.1) 


where 50 Q is the resistance of another leg of the 
Wheatstone bridge. Once the current is known, the power 
output of the sensor can be calculated using the following 
relation : 

P s = 1^ Rs (3.2) 

In order to negate the power supplied to the sensor 
due to heat losses to the substrate by conduction, a 
reference power level is calculated. The reference power 
level is calculated from the total bridge output voltage 
before the shock passage at time t = 0 . After the shock 

passes across the sensor, the result is purely the power 
supplied to the sensor due to changes in convective heat 

transfer across the sensor. Equation 3.3 illustrates the 
resulting power equation: 

APconv = Ps ~ P r ef (3-3) 


By neglecting changes radiation heat transfer during the 
experimental test time and negating the effects of power 
loss by conduction to the substrate, the resulting power 
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change, AP conv , is proportional to the convective heat 
transfer rate given by: 

q = hA(T w -T r ) (3.4) 


A best fit curve is placed through the transformed data 
using a computer software package?. Finally, t; he 

transformed data is reduced into a nond i mens 1 oria 1 

convective heat transfer parameter (i.e. Stanton number). 

The laboratory frame of reference is the physical 

coordinate system in which the experimental data is 

recorded and ultimately analyzed. Consequently, the 

Stanton number used in the experimental analysis is derived 

in the coordinate system fixed with respect to the 

laboratory. Thus, by applying the proportionality: 


APconv K qconv 


(3.5) 


and the definition of Stanton number, the resulting 
equation for convective heat transfer comparison becomes: 


St 


AP con v 

A T w — T rj^w^ u w ~ u e j ^ w 


(3.6) 


The Reynolds number (Eq. 2.20) is derived with respect to 
the laboratory reference frame and can therefore be used to 
produce a comparison of Stanton number as a function of 
Reynolds number. This becomes the most useful form in 
comparing the experimental data to theory and experimental 
work produced by Davies and Bernstein [6] and Roberts, et 
&L [13] . 
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The following chapter contains the oxperimetal 
results. These results are transformed into nond i mens i ona 1 
quantities for comparison with theory as we J 1 as for 
comparison with other experimental works [6,13] . Also 
obtained is a calculated heat transfer coefficient (h*) for 
.Judge’s [9] hot-film sensor computer simulation model. The 
derivation of the theoretical heat transfer coefficient is 
discussed in detail. Finally, the shocked gas pressure 
ratio data is discussed in the light of the effects of the 
diffuser downstream of the diaphragm. 


CHAPTER 4 


PRESENTATION OP RESULTS 
4.1 Experimental Results 

The data is acquired in terms of the voltage signals 
for the anemometer and the pressure transducer located at 
position 1 of the shock tube (see Fig. 3.2). Figure 4.1 
indicates the signal response of the hot-film sensor and 
position 1 pressure transducer at a driver pressure of 26.5 
psig. As stated in Section 3.2, the pressure signals at 

positions 1 and 2 are recorded to determine time of f 1 ight 
measurements of the passing shock wave. Furthermore, by 

recording the pressure response at position 1 , the 

experimentally measured pressure behind the normal shock 
wave can be used to deduce a shock wave speed to compare 
with the value found from the time of flight measurements 
[8] . Normal shock theory states that the pressure of the 
free stream gas behind the wave remains constant. From 
Fig. 4.1 (a) , it is evident that the pressure behind the 

shock wave decreases at the low driver pressures in the 
range of 37 to 41 percent from the time the normal shock 
passes across the hot-film sensor to the passage of the 
reflected shock (a period of approximately 3 ms). This 
percentage increases as the driver pressure is increased. 


36 







38 


At: the higher tested driver pressures, the decrease in 

pressure between the initial shock and the reflected shock 
varies from 47 to 51 percent. The theory used to correlate 
the experimental data is based on a constant free stream 
velocity behind the shock wave. Since the pressure 

response is showing a decrease from 37 to 51 percent over 
the test interval, this is an indication that the gas 
velocity is also decreasing over the test interval. The 
credibility of the theory is improved by identifying the 
experiment time interval as the first millisecond after 
shock passage across the hot-film sensor. 

Figure 4.2 illustrates a typical anemometer response 
at a low driver pressure and at a higher driver pressure. 
The response of the hot-film sensor is nearly horizontal 
and is characterized by a higher frequency voltage 
fluctuation at the higher driver pressures (Fig. 4.2 (b)) . 

The response of the anemometer at the lower driver 
pressures (Fig. 4.2 (a)) exhibits a response which 

quantitatively correlates with transient, compressible, 
turbulent boundary- 1 ayer theory over approximately the 
first millisecond of testing; there is otherwise little 
evidence of the expected power law development of a 
boundary layer behind the shock wave. This result seems to 
indicate that transition to turbulent flow is occurring 
very rapidly. The differences in response between the 
lower and higher driver pressures could be explained by 
less free stream t’rbulence occurring at the lower driver 
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Figure 4.2. Typical Anemometer R«ap° n »® at Driver Pressures 
9 of (a) 26.5 psig, and (b) 78.8 psig 
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pressures. At, the higher driver pressures, the boundary 
layer may be dominated by free stream turbulence which 
completely overcomes any ordered boundary-layer 

development . 

Because the pressure response is decreasing over the 


test t i me , the 

use of the first 

mil] 

i second more 

close 1 v 

approximates a 

constant pressure 

and 

constant free 

st ream 

velocity flow. 

Consequent 1 y , 

the 

co rre 1 at i on 

be t wet? n 


experiment and theory is limited to using the lower driver 
pressure tests over approximately the first millisecond of 
anemometer response. The leading negative impulse response 
seen in Fig. 4.2 is discussed in Section 4.2.2. 

The data set utilized for comparison with Davies and 
Bernstein [6] is narrowed to a set of six experimental 
tests. The data set consists of experiments 1-3, 22, 3 0, 
and 47. These tests were made on a 500 /is/div time base at 
the low range of driver pressures. The remaining tests are 
recorded on the smaller time base (20 /is/div) and/or at the 
larger driver pressures. Experiments 1-3 exhibit the same 
trends as the remaining three tests but are omitted from 
comparison with boundary- 1 ayer theory because the sensor is 
skewed approximately 60 degrees clockwise from the standard 
operating direction (see Fig. 3.4 for standard sensor 
orientation). 1 h i s causes a gradual ’’roll over” of the 
signal response for approximately the first 0.2 ms instead 
o f the initial peak voltage rise followed by a gradual 
decrease in volta b as seen in Fig. 4.2 (a). At the higher 



driver pressures. 


the ’’roll over” 
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at shock transit is more 
pronounced with the sensor in the «0 degree skewed 

orientation than at the lower driver pressures. Figure 4.3 
illustrates the effect of sensor orientation for a typical 
set of experimental data at the mid-range of tested driver 
pressures. The effect of sensor orientation is further 
tested by rotating the sensor in 90 degree increments from 
the normal flow orientation and monitoring the anemometer 
signal response. The changes observed in sensor response 

due to orientation are qualitative in nature. The sensor 
exhibits a ’’roll over” in anemometer response at the 

initial test time as discussed previously and an initial 
lag in anemometer response occurs as compared to the 
standard sensor orientation. This lag may be attributed to 
the misalignment of the hot-film sensor with respect to the 
pressure transducer location. The shock wave passes across 
the pressure transducer at position 1 before passing across 
the hot-film sensor. This leaves a data set for comparison 
with Davies [6] consisting of experiments 22, 39, and 47. 

' Figure 4.4 shows the anemometer responses obtained for each 
of these three experiments. Table 4.1 lists the relevant 
measured and calculated values of this data set. The shock 
velocity, u w is determined from time of flight 
measurements. The measured shock velocity is within 1.0 
percent of the calculated shock velocity using normal shock 
theory [8]. The ratio of the driven pressure to the 

undisturbed gas pressure (P 2 /P t ) is measured using the peak 
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Figure 4.3. Effect of Sensor Orientation on the Anemometer 
Response for, (a) Normal Operating Position, and 
(b) 60 Degree Skewed Orientation 










Table 4.1. Measured and Calculated Experimental Parameters 
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voltage signal obtained at the time the shock passes across 

the sensor. The calibration factor used for the pressure 

1 

transducers is 50 mv/psi. The function h*t 5 is the 

calculated theoretical convective heat transfer coefficient 
for a turbulent boundary layer and will be discussed in the 

following section. 

4.2 Heat Transfer Comparison 

4.2.1. Transient Stanton Number 

With the data set defined in Table 4.1, the Stanton 
number as a function of Reynolds number for each experiment 
in the data set is determined. The experimental Stanton 
number (Eq. 3.6) is calculated using the measured 
temperatures and the deduced velocities recorded at the 
time of the test. The remainder of the flow variables are 
calculated using normal shock theory. Also, the thermal 
and transport properties are calculated with respect to the 
operating temperature of the hot-film sensor since this 
temperature is taken to be the wall temperature (T w ) • The 
Reynolds number is calculated using Eq . (2.20) with the 

further modification of Eq . (2.21) for comparison with 

Dav ies [6] . 

The turbulent theory for compressible, transient 
boundary- layer flow suggests that the Stanton number is a 
function of Reynolds number to the -0.2 power (see Eq . 
(2.30)). Consequently, by placing the data on log-log 
scales and taking a least squares fit of order 1, a 


comparison can 


be made between theory and experiment. 


The 
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data is placed on log-log scales and a regression of order 
1 is made to obtain the coefficient and power of the best 
fit curve. A regression for experiment 47 is shown in Fig. 
4.5. Figure 4.5 (a) shows a regression for the range of 
data from the time the normal shock passes across the hot- 


film sensor to the 

return of the reflected s 

h ock 

wave wh i ch 

is approximately 

3 . 3 ms . 

F i gu re 

4 .5 

(»>> 

shows the 

regression for a 

segment 

of the dat 

a for 

1 ess 

than the 

first millisecond 

wh i ch 

cor re 1 ates 

with 

the 

tu rbu lent 

boundary- layer theory. 

Table 4.2 

lists 

the 

regress i on 

coefficients for 

a powe r 

1 aw fit ; 

3 l 1 ong 

with 

the t i me 


interval in which the regression is performed for three of 
the low pressure ratio experiments. Figure 4.6 shows the 
corresponding segment of experimental data and regression 
fits for the experiment numbers listed in Table 4.2. 

Table 4.2. Coefficients for Power Law Curve Fit* 


Exp . 
No. 

p</Pi 

Ao 

A, 

r 2 

T i me 
f rom 

(ms) 

to 

22 

2.82 

0.317 

-0. 171 

0.827 

0 . 0800 

0.800 

39 

2.79 

0.402 

-0.247 

0.879 

0.0700 

0.820 

47 

2.78 

0.279 

-0 . 183 

0.889 

0 . 0800 

0.760 


’Note: St = A 0 Rexs 


Table 4.2 indicates a satisfactory power law curve fit over 


the indicated time interval. 


The time 


interval for the 


47 
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Figure 4.5. 


Experimental Heat Transfer Correlation at a Driver 
Pressure of 26.5 psig for (a) All Data within Test 
Time, and (b) Data within First Millisecond 
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Figure 4.6. Stanton Number Power Low Regressions within the First 
Millisecond After Shock Transit 
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regression fit, begins with the peak voltage value recorded 
by the oscilloscope until the beginning of the first large 
intermittent voltage fluctuation which varies from 0.760 to 
0.820 ms after shock transit across the sensor. The 
correlation coefficent, r 2 , is in the 80 percentile range, 
which lends evidence for the deduction that a turbulent 
boundary layer is developing after the shock passage. 

The next step in the process of data comparison is to 
determine the theoretical Stanton number and compare this 
result with the experimental Stanton number. Using the 
equations developed in Section 2.4, the turbulent Stanton 
number as a function of Reynolds number is calculated. The 
theoretical transient Stanton number for experiment 47 is 
plotted along with the experimental Stanton number in Fig. 
4.7. It is seen from Fig. 4.7 that the experimental 
Stanton number is higher than theory predicts for the same 
Reynolds number. This trend is repeated for the other two 
test cases, which are not shown. Plotted with the 
theoretical and experimental heat transfer correlations in 
Fig. 4.7 is the Stanton number correlation which fits 
Davies [6] experimental data set. The experimental work by 
Davies [6] also showed evidence of a 1 am i nar- to- tu rbu 1 ent 
transition, an effect not seen in the eurren t work. Figure 
4.7 gives a qualitative comparison between the work of this 
thesis and experimental work produced by other researchers. 

A qualitative comparison is made in comparing the 
results of this work to the experiments conducted by 
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Roberts, £_t aj_ , [13] who used very thin thermocouples 
mounted in the shock tube wall. Roberts, rt aj_ showed heat 
transfer results as a function of time for test times of 
200 gs with good agreement with Mirels’ turbulent boundary- 
layer theory [4] . The shock velocities range from 1600 to 
1700 m/s which indicate that the driver pressures are much 
higher than those tested in this thesis. The shock 

velocities at the driver pressures tested in this thesis 
range from 382 m/s at the lowest driver pressures to 423 
m/s at the highest driver pressures. The ’’clean” data in 
[13] shows the same randomness and type of fluctuations as 
the data obtained in this thesis. Roberts, e£ [13] 

’’dusty” data exhibits chaotic behavior over the short test 
period. Consequently, the conclusion may be reached that 

flow obstructions, such as dirt or wall roughness in the 
shock tube, causes the signal response of the hot-film 
sensor to become chaotic and prevents good agreement with 
transient, compressible, turbulent boundary- 1 ayer theory. 

4.2.2. Heat Transfer Cpeff 

A transient convective heat transfer coefficient is 
calculated for use as the forcing function in Judge's 
computer simulation of a hot-film sensor mounted on an 
insulating foam substrate [9]. The computer simulation can 
then predict the anemometer response at a given set of flow 
conditions. This result may serve to calibrate the hot- 
fi lm sensor mod el. The t ran s i e n t he at transfer coefficient 


is defined by Eq . (2.19) for a laminar boundary layer. 
From analyzing the experimental data in Table 4.2, it is 
apparent that the data correlates with turbulent boundary- 
layer theory. Consequently, the turbulent convective heat 
transfer coefficient is utilized in the computer simulation 
in order to compare the simulated heat flux with the 
experimental result. Thus, Eq . (2.32) is used to determine 
the theoretical heat transfer coefficient. The theoretical 
heat transfer coefficient used in the computer simulation 
requires a modification to Eq . (2.32) . 

The temperatures that are experimentally recorded 
during the shock tube tests are the ambient air 
temperature, T^ , and the operating temperature of the hot- 
film sensor, T w . The convective heat transfer coefficient 
is calculated from the heat flux, where the temperature 
gradient is taken with respect to the recovery temperature, 
T r (see Eq . (2.16)). To maintain consistency between 
theory, experimental data, and the computer simulation, the 
temperature gradient for all three is taken with respect to 
the free stream temperature, T e . By adding and subtracting 
T w to Eq . (2.14); and, with further algebraic manipulation, 
the equation for the heat transfer coefficient ( Eq . (2.32)) 
takes the form: 


h* 


I 5 


0 . 0460 * kfn^m ( u w u e)~ 

r — 

p- 

L 





(0)Ec“ 

2 


(4.1) 
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Utilizing Eq . (4.1) as the convective heat transfer 

coefficient, consistency is maintained between the 

calculated convective heat transfer coefficient and the 
computer simulation model. 

As noted earlier, an interesting development in the 
determination of the heat transfer coefficient arises in 
approximately the first 25 /is of the anemometer response. 
A very pronounced dip in anemometer response occurs in this 
time range as seen by viewing the experimental data in Fig. 
4.8, which is recorded on a 20 microsecond time-base. The 
response of the hot-film sensor decreases to a negative 

peak voltage value and then increases sharply to a positive 
peak voltage value within the time frame of 25 /is . This 

drop in voltage increases as the driver pressure is 
increased. Figure 4.8 illustrates this increase in 

magnitude as the driver pressure is increased. Figure 3.3 
shows the sensor dimensions and the value of the 
characteristic length (0.0762 cm). The calculated time for 
the shock to pass across the sensor for Exp. 38 is 2.0 /is. 
In the DISA paper [5] , the shock front curvature was 
investigated by placing a hot-wire anemometer at varying 
distances from the shock tube wall. The results of this 
paper show a slight dip in the initial response of the hot- 
wire anemometer at a shock Mach number equal to 1.6. 
This is a qualitative verification that the voltage drop is 
characteristic of the anemometer response during passage of 



volts volts 
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a shock wave, but no works were found which clearly discuss 
this phenomenon. 

A hypothesis is generated in order to explain the 
phenomenon of this drop in voltage, which corresponds to a 
decrease in the convective heat transfer from the hot-film 
sensor. It is hypothesized that as the shock passes across 
the hot-film sensor, the temperature of the free stream gas 
behind the wave rises abruptly. This assumption is based 
upon normal shock theory. It is further hypothesized that 
the step change in free stream temperature across the 
shock wave precedes the development of the velocity 
boundary layer because of a mass transport lag. 
Consequently, natural convective conditions prevail prior 
to the build-up of the shocked gas velocity boundary layer, 
but the sensor-gas temperature difference is less. Thus, a 
smaller temperature gradient along with the same natural 
convective heat transfer coefficient produces a sudden 
decrease in the heat flux from the sensor. The current 

supplied to the hot-film sensor to maintain its constant 
operating temperature is decreased, which produces the 
initial drop seen in the experimental results for the first 
few microseconds of the test. At the start of the velocity 
boundary- layer growth the heat flux from the sensor is a 
maximum since the convective heat transfer coefficient is 
very large. As the velocity boundary layer grows, the heat 
flux from the sensor decreases because of the decrease in 
heat escaping the boundary layer, and the anemometer 
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supplies less current, to maintain the constant operat, jng 
temperature of the hot-film sensor. 

The transient convective heat transfer coefficient is 
divided into three parts. The first is an estimate of the 
amount of natural convection occurring before the start of 
the test from the initial power being supplied to the 
sensor to maintain its operating temperature. Second, a 
logical way of estimating the heat transfer coefficient 
from time t=0 to approximately t = 25/is for input into the 
computer simulation is made based on the prior hypothesis. 
Finally, the transient response derived from M i re 1 s [3] is 
used for times greater than measured shock transit times. 

The estimate of a natural convective heat transfer 

coefficient is made by assuming a characteristic length 

equal to the area of the sensor divided by the length of 

the sensor (0.0762 cm) . It is also assumed that this 

characteristic length can be associated with a small wire 

(horizontal cylinder) . With these assumptions, Morgan [14] 

defines a natural convective correlation for various ranges 

of Grashof numbers. The Grashof number is a dimensionless 

* 

parameter describing the ratio of buoyancy forces to 
viscous forces. After calcuating the Grashof number 

(approximately 0.12) and applying the correlation given by 
Morgan [14] , an estimate of the natural convective heat 
transfer coefficient is made yielding a value of 8.0 x 10 1 
W/m 2 -K, ± 25 percent due to the uncertainty in the 


cor re 1 at i on . 
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During the t i me that the out-put, voltage drops and 
r ises to its maximum value (0<t<25gs) , the heat transfer 
coefficient is assumed to fc>e a constant. The calculated 
time, using the deduced shock velocity, for the shock to 

pass across the sensor ranges from 1.8 to 2.0 ps at the 

tested driver pressures. The measured time half-way across 
the pressure pulse is on the order of 5 to 12 A* s • ^ 

discrepency exists between the calculated shock transit 
time across the sensor and the measured transit time. 
This discrepency may be attributed in part to the lag 
associated with the mass transport of the gas behind the 
shock. Another possibility relates to the response time of 

the anemometer. The transit time of the shock moving 

across the sensor is less by a factor of 10 than the 

frequency response of the anemometer (~ 100 kHz) . 

Consequently, the anemometer is not responding fast enough 
to capture this dip due to shock transit. 

Increasing the temperature from ambient conditions to 
the shocked gas temperature, and evaluating the heat 
transfer coefficient with respect to new flow properties, 

does not decrease the heat transfer coefficient 

significantly. Thus, the convective heat transfer 

coefficient is maintained at its natural convective value 
since the decrease is assumed to be within the error of 

previous assumptions. The transient portion of the 

convective heat transfer coefficient is derived from Eq. 

(4.1) for a turbulent boundary layer. By combining the 
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three portions, a. complete heat transfer coefficient is 
calculated for a given set of test conditions. Figure 4.9 
shows a convective heat transfer coefficient curve for Exp. 
47. Thus, the computer model can simulate the shock tube 
test cases using the combined heat transfer coefficient to 
predict a theoretical anemometer response . Tfi i s can 
ultimately be compared to the experimental anemometer 
response . 


4.3 Effects of Divergent Diffuser 
It is seen from Fig. 3.2 that the NASA Langley Shock 
Tube has a divergent diffuser located downstream of the 
diaphragm. Because of the diffuser, constant area shock 
tube theory cannot be used to predict the pressure (P 2 ) 
behind the wave as the wave travels down the shock tube. 
Consequently, modifications are made in the constant area 
shock tube theory to account for the effects of the 
diffuser. Bill Chapin, working for NASA Langley IRD, 
developed a correction in the existing shock tube theory 
using steady and unsteady isentropic relations along with 
applying the Rank i ne-Hugon i ot relations to account for the 
diffuser effect. Through a personal communication with Mr. 
Chapin, the correction to one-dimensional shock tube theory 
was obtained and used to estimate shocked gas and driver 
pressure ratios. Figure 4.10 shows the experimetal data 
for the series of shock tube tests. Also plotted is the 
correction in the constant area shock tube theory, which 
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ure 4.9. Transient Convective Heat Transfer Coefficient 



Constant Area Shock 



Figure 4.10. Shocked Gas Pressure Ratio Versus Driver Gas Pressure Ratio 
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accounts for the divergent diffuser, as well as the 
constant area theoretical curve. Fairly good agreement 
exists between the proposed correction and the test data. 

The following chapter contains the conclusions drawn 
from the results obtained from the shock tube tests and 
analysis. Recommendations are made which would extend the 
understanding of the transient response of hot-film 
anemometers • 



CHAPTER 5 


CONCLUSIONS AND RECOMMENDATIONS 

The transient response of the hot-fim anemometer has 
been analyzed with a combination of results. First* 
repeatable data is easily obtained at each of the tested 
driver pressure ranges. The experimental comparison of the 
data with transient turbulent compressible boundary- layer- 
theory seems to be valid for test times ranging from 0.2 to 
0.8 ms at driver pressures of approximately 26.0 psig. 
Even at these driver pressures with relatively good curve 
fits, the predicted Stanton number for turbulent flow is 
less than the Stanton number obtained from the experimental 
data. The slopes are similar (see Fig. 4.4), but the 
coefficients have different values. The experimental data 
correlates with turbulent boundary-layer theory within the 
first millisecond of testing. After the first millisecond, 
the anemometer response becomes intermittent. This may be 
due to flow disturbances occurring inside the shock tube 
due to rough walls, micro-cavities from the ports in the 
wal 1 s , or unsteady turbulent free stream flow. 

The comparisons with other experimental works [6,13] 
give reasonable qualitative comparisons. The analysis of 
the anemometer response at the higher driver pressures must 
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l y coni ro 1 1 ed 


Mounting the hot-film sensor on a surface in the center of 
the shock tube [6] might produce results that better 
correlate with theory over a wider range of driver 

pressures, since this will help minimize the flow 

disturbance due to the walls. Also, an anemometer with a 

higher frequency response may produce results which better 
correlate with compressible, boundary- 1 aye r theory. A 

surprising result of this work is the occurrence of reduced 
heat transfer as the shock wave passes the sensor. This 

effect has not been clearly demonstrated in the other 


experimental works. Further investigation is needed to 

fully understand this phenomenon. Only then will a 

reasonable means be available to predict the magnitude of 


this drop in voltage. 

The computer simulation of the anemometer response [9] 
is to be investigated in an extension of this project with 
hope that agreement between these experimental results and 

wi 1 1 tend to val idate the 


predicted anemometer response 
hot-film sensor model. 
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APPENDIX A 
LIST OF EQUIPMENT 
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Instrument 

1- P.C.B. Model 484B Line Power Unit 

2. P.C.B. Model 484B Line Power Unit 

3. Wallace and Tierman Pressure Gauge 
(150 ps i (1034 kPa) maximum) 

4. Fluke 2190A Digital Thermometer 

5. Fluke 2190A Digital Thermometer 

6. Mylar sheets for shock tube diaphragms 

7. NASA Langley Shock Tube 

8. Gould Digital Storage Oscilloscope, 
Type 4035 


NASA I . D . 
501702 
501703 

C44598 

431062 

427579 


138584 


9, Hewlett-Packard Digital X-Y Plotter, 
Model 7470A 


403822 


10. Hewlett-Packard X-Y Recorder, 

Model 7046A 532098 

11. Bendix Model 790 Microbarograph C01791 

12. Two P.C.B. Piezotronic High-Resolution 
Pressure Transducers, Models 112A21 
(sens it i v i ty = 50mv/ps i ) 

13. Anemometer 


a . 

Dantec Type 

55M01 Main Unit 

054789 

b. 

Dantec Type 

55M10 CTA Standard Bridge 


c . 

Dantec Type 

55M05 Power Pack 

054790 

d. 

5 ft. cable (connecting hot-film sensor 
to CTA standard bridge) 



Hot-film sensor mounted on an 
insulating foam substrate 



Instrument, Old Dominion 

14. Tektronix Digital Oscilloscope, 

Model 2430A 
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Un i versi ty I . D . 

32105 


15. Maxar 386 AT Computer with HP- IB Board 


33157 



APPENDIX B 

DATA ACQUISITION PROGRAM LISTING 


10 ’ Set up program for MS-DOS HP- IB I/O Library 

20 ’ For use independent of the PC instrument bus system 

30 DEF SEG 

40 CLEAR , &HFE00 

50 I=&HFE00 

60 ’PCIB. DIRS represents the directory where the library 
files are located 

70 PC IB .DIRS = ENVIRONS (” PCIB” ) 

80 1$ = PCIB. DIRS + ” \PCIBILC . BLD” 

90 BLOAD IS , &HFE00 

100 CALL I (PCIB. DIRS, 17., J7.) 

110 PCIB. SEG = 17. 

120 IF J7.=0 THEN GOTO 170 
130 PRINT ” Unable to load.” ; 

140 PRINT ” (Error #”;J7.;”)” 

150 STOP 

160 ’ Define entry points for setup routines 

170 DEF SEG = PCIB. SEG 

180 O.S =5 

190 C.S = 10 

200 I.V =15 

210 I .C =20 

220 L.P =25 

230 LD.FILE = 30 

240 GET. MEM = 35 

250 L.S =40 

260 PANELS = 45 

270 ’ Establish error variables and ON ERROR branching 

280 DEF. ERR = 50 

290 PCIB. ERRS = STRINGS (64 ,32) 

300 PCIB. NAMES = STRINGS ( 1 6 , 32) 

310 CALL 

DEF . ERR(PCIB . ERR, PCIB . ERRS , PCIB . NAMES , PCIB . GLBERR) 

320 PCIB . BASERR = 255 
330 ON ERROR GOTO 600 
340 J=-l 

350 I$=PCIB . DIRS+” \HPIB . SYN” 

360 CALL O.S (IS) 

370 IF PCIB . ERRoO THEN ERROR PC IB. BASERR 

380 ’ Determine entry points for HP-IB Library routines 

390 1=0 

400 CALL I . V( I , IOABORT, IOCLEAR, I0C0NTR0L, IOENTER) 

410 IF PC IB. ERRoO THEN ERROR PC IB. BASERR 
420 CALL I .V(I , IOf NTERA, IOENTERS, I0E0I , I0E0L) 
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IF PCIB.ERRoO THEN ERROR PCIB.BASERR 
CALL I . V ( I , IOGETTERM, IOLLOCKOUT , IOLOCAL, IOMATCH) 

IF PCIB.ERRoO THEN ERROR PCIB.BASERR 

CALL I . V ( I , IOOUTPUT, IOOUTPUTA, IOOUTPUTS, IOPPOLL) 

IF PCIB.ERJloO THEN ERROR PCIB.BASERR 
CALL I . V ( I , IOPPOLLC , IOPPOLLU, IOREMOTE , IORESET) 

IF PCIB.ERRoO THEN ERROR PCIB.BASERR 
CALL I . V ( I , I OSEND , I OSPOLL , I OSTATUS , I OT I MEOUT ) 

IF PCIB.ERRoO THEN ERROR PCIB.BASERR 
CALL I . V ( I , I OTR I GGER , I ODMA , J , J ) 

IF PCIB.ERRoO THEN ERROR PCIB.BASERR 
CALL C.S 

I $=PC I B . D I R$+” \HP I B . PLD” 

CALL L.P(I$) 

IF PCIB.ERRoO THEN ERROR PCIB.BASERR 
GOTO 680 

> Error handling routine 
IF ERR=PC I B . BASERR THEN GOTO 630 
PRINT ’’BASIC error #”;ERR;” occurred 
STOP 

TMPERR = PC I B . ERR 

IF TMPERR = 0 THEN TMPERR = PCIB.GLBERR 
PRINT ”PC Instrument error #”;TMPERR;” detected at 

M . 

PRINT ’’Error: ”;PCIB.ERR$ 

STOP 

COMMON PC I B . D I R$ , PC I B . SEG 
COMMON LD . F I LE , GET . MEM , PANELS , DEF . ERR 

PCI B . BASERR , PC I B . ERR , PC I B . ERRS , PC IB. NAMES, PC IB. GLBERR 
Tin pnMMON 

IOABORT IOCLEAR I0C0NTR0L, IOENTER, IOENTERA , IOENTERS, I0E0I , 10 
EOL IOGETTERM IOLLOCKOUT , ioLOCAL, IOMATCH , IOOUTPUT , IOOUTPUTA . 
I OOUTPUTS , I OPPOLL , I OPPOLLC , I OPPOLLU , IOREMOTE , I ORESET , I OSEND , 
I OSPOLL , I OSTATUS , I OT I MEOUT , I OTR I GGER , I ODMA 


430 
440 
450 
460 
470 
480 
490 
500 
510 
520 
530 
540 
550 
560 
570 
580 
590 
600 
610 
620 
630 
640 
650 
1 i ne 
660 
670 
680 
690 
700 


l n 


1 i ne ” ; ERL 


720 

730 

740 

750 

760 

770 

780 

790 

800 

810 

820 

830 


FALSE 
TRUE 
NOERR 
EUNKNOWN 
ESEL 
ERANGE 
ETIME 
ECTRL 
EPASS 
ENUM 
EADDR 


= 0 


NOT FALSE 
0 

100001 ! 

100002 ! 

100003! 

100004! 

100005! 

100006! 

100007! 

_ COMMON FALSE°°TRUE, NOERR, EUNKNOWN, ESEL, ERANGE, 
ETIME, ECTRL, EPASS, ENUM, EADDR 
840 ’ End Program Set-up 

850 ? User program begins , . 

860 INPUT ’’Enter the file name with path for storing 

”, RESFILES 

870 ’ Open the file to download the data from scop 



880 OPEN RESFILES FOR OUTPUT AS #1 

890 INPUT ’’Enter the EXPERIMENT number : ”, INFOS 

900 PRINT #1, TAB (5) "Experiment # : INFOS 

910 PRINT ’’User Program For TEKTRONIX 2430A Digital 
Oscilloscope ” 

920 ’ Set the addresses for GPIB card, Scope and Plotter 

930 ISC = 7 : SCOPE = 13 : PLOTTER = 5 

940 CRT = ISC * 100 + SCOPE 

950 PLT = ISC * 100 + PLOTTER 

960 ’ Initialize the variables 

970 MAX. ELEMENTS = 1024 : ACTUAL . ELEMENTS = 0 

980 MAX. LENGTH = 255 : ACTUAL . LENGTH = 0 

990 OPTION BASE 1 

1000 ’ Define, the Array dimensions 

1010 DIM WFM ( 1 024 ) 

1020 DIM VOLTAGE (1050) 

1030 ’ Reset the GPIB card by calling IORESET command 

1040 CALL IORESET (ISC) 

1050 TIMEOUT = 5 

1060 ’ Clear the Scope and set a time for the device to 

respond in secs. 

1070 CALL I OCLE A R ( CRT ) 

1080 CALL I OT I MEOUT ( I SC , TIMEOUT) 

1090 ENABLE = 1 

1100 CALL I0E0 I (ISC, ENABLE) 

1110 ’ Set the GPIB to remote 
1120 CALL I OREMOTE ( I SC) 

1130 CALL IOCLEAR(ISC) 

1140 ’TRIGGERING COMMANDS basically for setting preset 
trigger points using 

1150 ’the Atrigger and Btrigger positions 

1160 C0MM$ = ’’atrigger sou rce : chi , pos i t i on : 04” : LENGTH = 

LEN (COMMS) 

1170 CALL I OOUTPUTS (CRT, COMMS, LENGTH) 

1180 COMMS = ’’btrigger sou rce : ch 1 , pos i t i on : 04” : LENGTH = 

LEN (COMMS) 

1190 CALL I OOUTPUTS (CRT, COMMS, LENGTH) 

1200 COMMS = ’’PATH OFF” : LENGTH = LEN (COMMS) 

1210 CALL I OOUTPUTS (CRT, COMMS, LENGTH) 

1220 ’Settings of the Scope 
1230 ’ VERTICAL Commands set up 

1240 COMMS = ” VMODE CH 1 : ON , CH2 : ON ; BWLIMIT FULL” : LENGTH 
LEN (COMMS) 

1250 CALL I OOUTPUTS (CRT, COMMS, LENGTH) 

1260 COMMS = ” CH 1 VOLTS : 100E-3 , VAR : 0 .COUPLI NG : AC ; CH2 
VOLTS : 200e -3, VAR :0, COUPLING: AC” : LENGTH = LEN (COMMS) 

1270 CALL I OOUTPUTS (CRT, COMMS, LENGTH) 

1280 ’ HORIZONTAL Commands set up 

1290 COMM$=” HOR ASECDI V : 20E-6 , BSECDI V : 20E- 

6, POSITION: 1 .28E+2” : LENGTH = LEN (COMMS) 

1300 CALL I OOUTPUTS (CRT, COMMS, LENGTH) 

1310 ’ WAVEFORM Commands set up 

1320 COMMS = ’’START l; STOP 1024; DATA ENCDG : ASCI I” : 
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LENGTH = 


LENGTH = LEN (COMMS) 

1330 CALL I OOUTPUTS (CRT, COMMS, LENGTH) 

1340 ’ Acquisition commands 

1350 COMMS = ’’RUN ACQUIRE; ACQUIRE MODE: NORMAL” 

LEN (COMMS) 

1360 CALL I OOUTPUTS (CRT, COMMS, LENGTH) 

^Start~Acqu i ri ng the DATA from SCOPE to COMPUTER 
IF (CHANL = 1) THEN PRINT ’’Reading CHANNEL 1 
IF (CHANL = 2) THEN PRINT ’’Reading CHANNEL 2 
PRINT #1, ’’Readings of Channel Number : CHANL 

IF CHANL = 1 THEN GOTO 1440 ELSE GOTO 1450 
’Read the values from channel 1 or 2 depending upon 
channel 


1370 

1380 

1390 

1400 

1410 

1420 

1430 

the 


1440 
1460 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
1620 
#1 , : 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
TFL 
1730 
1740 
1750 
1760 


COMMS = ’’DATA SOURCE: CHI” : LENGTH = LEN (COMMS) 

LEN (COMMS) 


GOTO 


COMMS = ’’DATA SOURCE :CH2” : LENGTH = 

CALL I OOUTPUTS (CRT, COMMS, LENGTH) 

’ Start getting the data values from scope 
COMMS = ’’CURVE?” : LENGTH = LEN (COMMS) 

CALL I OOUTPUTS (CRT, COMMS, LENGTH) 

CALL I OENTERA ( CRT , WFM ( 1 ) , MAX . ELEMENTS , ACTUAL . ELEMENTS) 
COMMS = ’’WFMPRE? YOFF” : LENGTH = LEN (COMMS) 

CALL I OOUTPUTS ( CRT , COMMS , LENGTH ) 

CALL IOENTER(CRT , YOFF) 

PRINT #1, ” Yof f read as : YOFF 

COMMS = ’’WFMPRE? YMULT” : LENGTH = LEN (COMMS) 

CALL I OOUTPUTS (CRT, COMMS, LENGTH) 

CALL IOENTER(CRT , YMULT) 

PRINT #1, ” Ymu 1 t read as : ”, YMULT 

COMMS = ’’WFMPRE? XI NCR” : LENGTH = LEN (COMMS) 

CALL I OOUTPUTS (CRT, COMMS, LENGTH) 

CALL IOENTER(CRT , XINCR) 


PRINT 
PRINT 
PRINT 
PRINT 
’ Get 
FOR I 


#1 , 
#1 , 
#1 , 
# 1 , 
the 
= 1 


'Xincreraent is read as 


’ Po i nt 


T i me 


, XINCR : PRINT 
Voltage” 


1024 data points read by the scope 
TO 1024 

VOLTAGE ( I ) = (WFM ( I ) - YOFF) * YMULT 
PRINT #1, I TAB(12) (I-1)*XINCR TAB(26) 

NEXT I 

CHANL = CHANL + 1 
IF CHANL < 3 THEN GOTO 1380 


VOLTAGE ( I ) 


IT Ulinnu n u - — - — 

INPUT ’’Enter the TIME of FLIGHT from GOULD Scope 


PRINT #1 
PRINT #1 
CLOSE #1 
’ To plot the 
the instructions 
1770 PRINT ” 

1780 PRINT : PRINT 


TIME of FLIGHT 
TEKTRONIX 2430 A 


: ” , TFL 
Scope Read i ngs’ 


curves on 


the plotter from scope follow 


NOTE 



To PLOT the waveforms follow the instructions 



1790 PRINT ” 


1800 PRINT ”1. Make sure NOT to touch the DIP switches of 
PLOTTER” 

1810 PRINT ”2. Insert a blank paper in the plotter” 

1820 PRINT ”3. Remove the GPIB cable from the COMPUTER” 

1830 PRINT ”4. Using the OUTPUT button make the SCOPE to 
DEVICES Mode using the SETUP Bezel” 

1840 PRINT ”5. Again press the OUTPUT button and press PLOT 
bezel” 

1850 CALL I 00UTPUTS( CRT, COMMS, LENGTH) 

1860 COMMS = ”?7.” : LENGTH = LEN(COMMS) 

1870 END 
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